最 小 妨 忌 公平 性 区 位 模型 及 案例 分 析 
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摘要 : 区 位 模型 广泛 应 用 于 公共 健康 、 义 务 教 育 、 应 急 管 理 和 物流 配送 等 领域 。 然 而 ， 主 流 区 位 问题 多 以 
服务 效率 为 目标 ， 注 重 设 施 成 本 、 距 离 成 本 、 履 盖 客 户 数量 等 指标 ， 未 考虑 服务 公平 性 。 部 分 区 位 模型 考 
虑 服务 空间 公平 性 ， 但 存在 公平 指标 多 且 选 择 使 用 缺乏 共识 、 公 平 指标 会 严重 扭曲 效率 指标 、 模 型 针对 简 
单 应 用 场景 和 缺乏 通用 性 等 局 限 。 本 文 创新 性 地 引入 空间 妒忌 这 一 概念 , 构建 公共 服务 最 小 嫉妒 目标 函数 ， 
并 将 其 整合 在 容量 约束 P 中 值 问题 (Capacitated P-Median Problem，CPMP) 模 型 目标 函数 中 ， 形 成 最 小 
妒忌 公平 性 区 位 模型 (CPMP-envy)。 新 模型 克服 了 最 小 方差 目标 会 扭曲 效率 指标 的 局 限 ,， 实现 公平 指标 与 
效率 指标 相 兼 容 。 本 文 考虑 城市 和 农村 地 区 人 口 空 间 分 布 特征 ， 以 及 数据 规模 等 因素 ， 使 用 三 个 典型 区 域 
数据 探索 经 典 区 位 模型 和 公平 性 区 位 模型 的 基本 特征 。 案 例 实验 表明 : 最 小 妒忌 目标 能 够 显著 地 改进 设施 
区 位 空间 公平 性 指标 ， 特 别 是 设施 数量 较 少 时 ， 公 平 性 指标 能 够 得 到 较 大 幅度 改进 。 最 小 妒忌 目标 在 公共 
服务 设施 布局 规划 中 具有 实用 性 ， 对 于 构建 公平 性 区 位 模型 具有 理论 价值 。 
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区 位 问题 是 地 理学 、 数 学 、 区 域 科学 、 管 理科 学 、 交 通 科 学 、 计 算 机 科学 等 学 科 共 同 关 
注 的 一 个 基础 性 问题 申 。 自 上 世纪 50 年 代 以 来 ， 区 位 问题 模型 建构 、 数 学 性 质 与 求解 算法 
研究 持续 不 断 。 经 典 的 区 位 问题 包括 区 位 集 窗 新 问题 (Location Set Covering problem, LSCP)?!, 
P 中 值 问题 (P-Median Problem, PMP)Pl, P 中 心间 题 (P-Center problem, PCPA, RAH ss KA 
问题 (Maximum Covering Location Problem, MCLP)P!, Wit X (iz [i] ii (Uncapacitated Facility 
Location Problem, UFLP)I 等 。 随 区 位 问题 研究 的 深入 和 应 用 范围 的 扩大 ， 经 典 问 题 中 逐渐 
加 入 了 设施 容量 限制 ， 如 容量 约束 P 中 值 问题 (Capacitated P-Median Problem, CPMP)", 4 
量 约束 P 中 心 问题 (Capacitated P-Center problem, CPCP)ISI, 725 R px fr E 78 si [Al 
(Capacitated Location Set Covering problem, CLSCP) 和 容量 约束 最 大 覆盖 区 位 问题 
(Capacitated Maximum Covering Location Problem, CMCLP)P!, 4 Œ A) IR v yi pc fiz qn] i 
(Capacitated Facilities Location Problem, CFLP)09 等 。 

上 述 问题 均 是 从 候选 区 位 中 选择 部 分 区 位 , 用 于 建设 服务 设施 , 并 确定 设施 为 哪些 需求 
提供 服务 。 其 中 ，PMP、PCP 和 MCLP 模型 均 需 事先 确定 设施 数量 ，PMP 最 小 化 总 距离 成 
AS, PCP 最 小 化 最 大 距离 ， 而 MCLP 最 大 限度 地 者 新 需求 ; LSCP 和 UFLP 不 需要 事先 确定 
设施 数量 , 前 者 尽 可 能 使 用 最 少 的 设施 或 者 总 成 本 最 低 的 设施 , 后 者 兼顾 设施 总 成 本 和 总 距 
离 成 本 。 上 述 模 型 广泛 应 用 于 公共 健康 、 义 务 教育 、 应 急 管 理 、 物流 配送 、 休闲 娱乐 等 领域 。 
然而 , 主流 区 位 问题 多 以 服务 效率 为 目标 , 注重 设施 成 本 、 上 距离 成 本 、 履 盖 客 户 数量 等 指标 ， 
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未 考虑 服务 公平 性 。 

公平 性 区 位 是 具有 挑战 性 的 一 类 区 位 问题 。 随 经 济 社会 快速 发 展 , 区 位 问题 已 广泛 应 用 
于 公共 服务 设施 规划 ， 且 越 来 越 注重 空间 公平 性 055。 国 家 《“ 十 四 五 ?公共 服务 规划 》 对 服 
务 设施 布局 提出 了 基本 要 求 : 需求 全 和 窗 新 ,质量 全 达标 ; 均等 剖 有 ， 便 利 可 及 。 为 满足 公共 
服务 设施 空间 公平 性 需求 ， 学 者 提出 多 个 公平 性 指标 ， 如 距离 极 值 、 方 差 、 变 异 系 数 、Gini 
系数 等 09， 并 基于 公平 性 指标 建立 各 种 公平 性 区 位 模型 。 

这 些 模型 可 归纳 为 两 大 类 。 第 一 ， 使 用 公平 性 目标 函数 。 将 经 典 区 位 模型 的 效率 目标 蔡 
换 为 公平 性 目标 函数 ， 如 最 小 化 距离 基尼 系数 0718、 最 小 化 B-mean 距离 (|、 最 小 化 距离 与 
期 望 值 偏差 中。Berman 使 用 最 小 化 距离 方差 目标 函数 ， 并 引入 平均 距离 约束 条 件 ， 防 止 公 
平 性 目标 扭曲 服务 效率 指标 PH。 王 法 辉 等 发 展 了 两 步 优化 法 (Two-step Optimization for Spatial 
Accessibility Improvement，2SO4SAT)， 先 使 用 MCLP 或 PMP 模型 选择 设施 区 位 ， 再 以 最 小 
化 2SFCA 可 达 性 方差 为 目标 分 配 资源 P229。 第 二 ， 使 用 多 目标 优化 模型 。 优 化 目标 包括 效 
率 指 标 函 数 和 公平 指标 函数 ,模型 求解 获得 Pareto 最 优 解 集 。 常见 的 方法 是 在 经 典 区 位 模型 
中 增加 最 小 化 距离 方差 、 最 小 化 覆盖 率 差 异 、 最 小 化 偏差 等 公平 性 目标 ,构造 多 目标 优化 模 
型 27394。 也 可 以 将 公平 指标 函数 与 效率 指标 函数 加 权 相 加 为 单 目 标 函 数 ， 例 如 ， 将 B-mean 
距离 和 平均 距离 加 权 之 和 作为 目标 函数 ， 使 新 冠 检测 点 选 址 兼顾 到 可 及 性 和 公平 性 093031, 

基于 公平 性 指标 所 建立 的 各 种 公平 性 区 位 模型 , 存在 较 大 的 局 限 性 , 严重 影响 其 应 用 成 
效 。 首 先 ， 公 平 性 指标 多 ， 相 互 差异 较 大 ， 如 何 选 择 指标 缺乏 共识 43。 其 次 ， 常 见 公平 性 目 
标 ( 如 方差 、 变 异 系数 和 Gini 系数 ) 会 严重 扭曲 效率 指标 089, 计算 结果 难以 获得 实用 的 设施 布 
局 方案 。 第 三 ， 公 平 性 指标 往往 是 非 线性 的 呈 ， 导 致 模型 构建 困难 ， 且 模型 计算 复杂 度 极 
高 。 第 四 ， 现 有 公平 性 区 位 模型 多 是 应 对 简单 场景 (如 未 考虑 设施 容量 )， 针 对 具体 应 用 ， 通 
用 性 不 足 。 例 如 ，2SO4SAI 方法 应 用 案例 较 多 ， 主 要 用 于 公平 地 分 配 资源 到 设施 区 位 ， 并 非 
选择 公平 性 区 位 。 

综 上 , 现 有 公平 性 区 位 模型 仍 存 在 较 大 的 局 限 性 ， 相 关 研 究 进 展 不 大 中 。 针 对 公平 性 区 
位 模型 面临 的 问题 与 挑战 ,本文 梳理 相关 研究 进展 , 综合 考虑 城市 和 农村 地 区 人 口 空间 分 布 
特征 ,以 及 数据 规模 等 因素 , 使 用 三 个 典型 区 域 数据 探索 经 典 区 位 模型 和 公平 性 区 位 模型 的 
基本 特征 。 通过 案例 实验 总 结 最 小 化 方差 区 位 模型 的 特征 及 局 限 , 并 尝试 构建 新 的 公平 性 区 
位 模型 .新 模型 试图 克服 现 有 模型 的 局 限 , 也 为 公平 导向 公共 服务 设施 布局 规划 提供 新 方法 。 


2 公平 性 区 位 模型 构建 


令 集 合 1={1,2.…nn} 表示 n 个 候选 服务 设施 区 位 ， 设 施 i(i e D 有 最 大 服务 容量 s, 和 国定 
成 本 fi; 集合 /三 {1,2… 区 表示 m 个 空间 单元 ， 单 元 j(j eJ) 的 服务 需求 量 为 wj， 空 间 单元 i 
与 单元 j 之 间 的 距离 为 dj 、 距 离 成 本 为 cy 。 定 义 布尔 型 决策 变量 xi 为 是 否 指派 设施 了 服务 
需求 单元 j; 布尔 型 决策 变量 y, 表 示 是 否 在 区 位 i 建设 施 。 基 于 这 些 定义 ， 可 以 构建 经 典 区 
位 模型 和 公平 性 区 位 模型 。 


CPMP 模型 : 
Minimize: Dicer Xie] W;dijXij (1 = 1) 
Subject to: Vier Xij = 1,Vj EJ (1-2) 
È jej WjXij < Siyi Vi EI (1-3) 
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dier Vi = P (1—4) 
Xip yi e {0,1}, vi E IjEJ (1-5) 
方差 最 小 区 位 模型 (Minimum Variance Location Problem, MVLP): 


Minimize: Ye; w;d; + Vie w;d'^ — 2X je; W;d;d' (2— 1) 
Subject to: Vier xij = 1, Vj EJ (2 一 2) 
2jeJwixXij < Siyo Vi € I (2 一 3) 

2ieryi = P (2-4) 

dj = Vier dijxij, Vj EJ (2 一 5) 

d'= Xjejw;dj / XMjejWj (2-6) 

Xip Yi €10,15 Vi E Lj EJ (2-7) 

dd > 0,vi EJ (2-8) 


MVLP 模型 中 ，dj 表 示 和 需求 点 j 使 用 服务 的 出 行距 离 ，d' 为 平均 出 行距 离 ， 目 标 (2-1) 为 出 行 
距离 方差 。MVLP 模型 是 一 个 二 次 规划 模型 ， 与 CPMP 模 相 比 ， 增 加 了 决策 变量 qj 和 d'， 也 
增加 了 相应 的 约束 条 件 。 

给 定 距离 均值 dq, 并 增加 最 大 均值 约束 ,可 将 MVLP 模型 转化 为 线性 规划 模型 (Minimum 
Variance Location Problem with Maximum Distance, MVLP-MD): 


Minimize: Xici 2 jej wi (di 一 d') xij (3— 1) 
Subject to: Vier xij = 1, Vj EJ (3 —2) 
YjejWwjxij < Siyo Vi € I (3—3) 
dier Vi = P (3-4) 

dier jez W;jdijxij € d' * Vie Wi (3 一 5) 
Xip Yi E {0,1}, Vi E Lj EJ (3 — 6) 


给 定 不 同 的 距离 均值 4 ，MVLP-MD 模型 能 够 获得 不 同 的 区 位 方案 。 

因 距 离 方 差 公平 性 指标 会 严重 扭曲 效率 指标 ， 本 文 将 最 小 方差 指标 修订 为 最 小 “ 妨 忌 ” 指 
标 。 给 定 距离 参数 d*， 假 定 茶 一 区 位 人 群 使 用 服务 出 行距 离 超 过 qd* 时 ,会 有 不 公平 感知 , 或 
者 会 对 使 用 服务 出 行距 离 小 于 d* 的 人 群 产 生 妨 忌 。 定 义 区 位 j 妒忌 值 为 : 


0, vd; < d" 
- 4—1 
: M - d), vdij > d' (D 
AP eo; RIZ" P Hb ERR: 
HEN «MA 
Minimize die WV; = Dier dies died" w;(dij -d ) Xij (4 = 2) 


目标 函数 (5-2) 仅 计算 出 行距 离 大 于 d* 人 群 的 距离 差异 平方 值 , 有 利于 降低 设施 布局 的 不 公平 
性 。 因 出 行距 离 较 远 的 需求 者 才 有 不 公平 感知 ， 该 指标 与 空间 可 及 性 指标 不 冲突 ,不 会 扭曲 
可 及 性 指标 。 

将 函数 (5-2) 加 入 CPMP 得 到 改进 的 区 位 模型 ， 记 为 CPMP-envy 模型 : 


. . . * 2 
Minimize: die dies Wjdijxij + Q diel È jeJ,dijzd* w;(dij -d ) Xij 5-1) 
Subject to: Vier xij = 1, Vj EJ (5 一 2) 
ies WjXij < Siyi Vl eI (5 —3) 
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Yiery; = P (5—4) 
Xij Yi € (015, Vi e Lj EJ (5 — 5) 
CPMP-envy 模型 中 ， 妨 忌 目 标 权 重 a 的 取 值 较为 重要 。 若 权重 a 数值 过 小 ， 妨 忌 目 标 不 
起 作用 , 模型 趋 近 于 CPMP 模型 ; 看 权重 w 数 值 过 大 , 妒忌 目标 会 扭曲 距离 成 本 目标 。 因 此 ， 
权重 a 数值 选择 的 原则 是 妒忌 目标 和 距离 成 本 目标 数值 接近 ， 不 会 扭曲 成 本 目标 ， 又 能 增 
加 空间 公平 性 。 
假设 距离 均值 为 d*, 距离 均 方 差 为 6 = yd*， 目 标 函 数 (5-D) 中 ， 第 一 部 分 数值 为 jeJ w; * 
qd*， 第 二 部 分 估 值 为 (0.5y?d”?) * Xjejw;. Sa = 2y ?/d', 目标 函数 (5-1) 中 ， 两 部 分 量 纲 一 
致 ， 且 数值 接近 。 注 意 ，y = 6/d* 是 距离 变异 系数 ， 当 y=0.5 I, æ = 8.0/d*， 当 y=0.6 时 ， 
a x 5.6/d*。 较 为 准确 地 估算 d* 和 y， 有 利于 CPMP-envy 模型 获得 较为 公平 的 公共 服务 。 


3 模型 实验 


3.1 实验 设计 


本 文 使 用 三 个 区 域 数据 探索 经 典 区 位 模型 和 公平 性 区 位 模型 的 基本 特征 ,三 个 区 域 特征 
WR 1 所 示 ， 考 虑 到 了 城市 和 农村 区 域 ， 也 考虑 到 了 数据 规模 。 区 域 GY 和 ZZ 将 人 口 数量 
作为 需求 量 ， 区 域 ZY 将 小 学 生 人 数 作为 需求 量 。 区 域 GY 选择 初级 中 学 、 乡 镇 居民 点 、 街 
道 办 事 处 ， 以 及 知 干 规模 较 大 行政 村 作为 候选 设施 区 位 ; 区 域 ZY 将 所 有 需求 点 作为 候选 设 
施 区 位 ; 将 区 域 ZZ 划分 为 0.7km2 正六 边 形 格 网 ,从 每 个 格 网 中 距离 中 心 位 置 最 近 的 6 个 需 
求 点 中 随机 选择 一 个 点 作为 候选 区 位 。 

表 1 案例 基本 特征 


Tab. 1 Summary of instance attributes 


区 域名 称 地 理 背 景 面积 (km?) 需求 点 数量 总 需求 数量 候选 设施 点 数量 设施 容量 (人 数 ) 


GY 农村 998.7 1276 819812 135 40000 
ZY 城市 13.4 324 3873 324 400 
ZZ 城市 297.8 6752 4223997 444 50000-90000 


所 有 区 域 案 例 中 , 距离 dij 采 用 干 米 单位 , 距离 成 本 crj = 1.0 * dijo CPMP 和 CPMP-envy 
模型 使 用 数学 启发 算法 求解 6B9。 对 于 每 一 个 案例 的 计算 结果 ,统计 常见 服务 指标 ， 如 设施 数 
量 、 平均 距 离 、 最 大 距离 、 方 差 、 标 准 差 、 变 异 系数 、Gini 系数 等 指标 。 令 dj = Vier dijxij: 
d = Dies Djey dijxij/ 了 jejwj， 距 离 方 差 (Van) 计 算 公式 见 (6-1)。 


Var = Xjepwidid) wildi- 可 (6 — 1) 
Yjejwj 
距离 标准 差 (Stdev) 和 变异 系数 (C 隐 见 公 式 (6-2) 和 (6-3): 

Stdev = Var ? (6 — 2) 
CV — Stdev/d (6 — 3) 

考虑 到 各 区 位 需求 量 ，Gini 系数 (GO 计算 公式 见 (6-3)。 

Yic] Lies Wiw j|di-dj 

ae m AEEA (6-9 
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3.2 案例 计算 结果 


案例 GY 计算 结果 统计 见 表 2。 表 中 各 栏目 依次 为 模型 名 称 、 模 型 参数 、 设 施 数量 及 各 
区 位 出 行距 离 的 均值 、 极 值 、 方 差 、 标 准 差 、 变 异 系数 和 Gini 系数 。 对 于 CPMP 模型 (P=23)， 
平均 距离 、 最 大 距离 、 方 差 、 标 准 差 、 变 异 系 数 和 Gini 系数 依次 为 1.948. 8.218. 2.318. 
1.522. 0.782 和 0.439; 增加 最 小 妒忌 目标 ， 平 均 距 离 增加 5.3%， 但 最 大 距离 、 方 差 、 标 准 
差 、 变 异 系数 和 Gini 系数 减 小 显著 ， 分 别 减 小 14.8%、21.0%、11.1%、15.7% 和 14.6%。 设 
施 数量 减少 (P=22) 或 者 设施 数量 增加 (P=24)， 除 平均 距离 指标 外 ，CPMP-envy 结果 的 最 大 距 
离 、 方 差 、 标 准 差 、 变 异 系数 和 标准 差 指 标 均 优 于 CPMP 模型 。 上 述 结果 表明 ， 增 加 最 小 妒 
鼠 目 标 ， 公 平 性 指标 可 明显 改善 。 

表 2 案例 GY 计算 结果 统计 


Tab. 2 Computational results from instance GY 


标准 差 变异 Gini 


模型 名 称 模型 参数 -- ee CE 
CPMP P-22 22 2.064 8.641 2.419 1.555 0.754 0.422 
CPMP-envy P=22,d*=2.0km,a=4.0 22 2.151 8.218 1.988 1.410 0.655 0.370 
CPMP P-23 23 1.948 8.218 2.318 1.522 0.782 0.439 
CPMP-envy P=23,d*=2.0km,a=4.0 23 2.052 7.002 1.831 1.353 0.659 0.375 
CPMP P=24 24 1.875 8.218 2.190 1.480 0.789 0.445 
CPMP-envy P=24,d*=2.0km,a=4.0 24 2.007 6.372 1.618 1.272 0.634 0.360 


fik: CPMP: 容量 约束 P 中 心 问 题 ， CPMP-envy: 最 小 妒忌 公平 性 区 位 模型 

1 为 CPMP 和 CPMP-envy 模型 (P=23) 计 算 结果 示意 图 。 从 统计 指标 看 ， 两 个 模型 结 
果 中 平均 出 行距 离 差异 不 算 太 大 (5.3%)， 后 者 方差 和 Gini 系数 指标 有 较 大 幅度 的 下 降 。 同 
时 ， 两 个 模型 所 选择 设施 区 位 差异 巨大 ，23 个 设施 区 位 中 ， 仅 11 个 位 置 相同 。 


1 GY 案例 CPMP( 左 ) 和 CPMP-envy( 右 ) 模 型 计算 结果 示意 
Fig.1 CPMP (left) and CPMP-envy (right) solutions from instance GY 


案例 ZY 计算 结果 统计 见 表 3。 表 中 各 栏目 含义 与 表 1 相同 。 对 于 CPMP 模型 (P=11)， 


FHEA, RAER, JA. PE, ARAM Gini 系数 依次 为 0.408、1.053、0.046、 
0.214. 0.526 和 0.285; 增加 最 小 妒忌 目标 ， 平 均 距 离 增加 2.9%， 但 最 大 距离 、 方 差 、 标 准 
差 、 变 异 系数 和 Gini 系数 减 小 显著 ， 分 别 减 少 4.596. 26.196. 14.096. 16.5991 17.3%。 对 于 
CPMP 模型 (P=12)， 增 加 最 小 妒忌 目标 ， 最 大 距离 、 方 差 、 标 准 差 、 变 异 系 数 和 Gini 系数 均 
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减 小 显著 。 
表 3 案例 ZY 计算 结果 统计 


Tab. 3 Computational results from instance ZY 


模型 模型 参数 umo Ra ya WEE 22 Qa 
CPMP P=11 11 0.408 1.053 0.046 0.214 0.526 0.295 
CPMP-envy P=11,d*=0.4km,a=20.0 11 0.420 1.006 0.034 0.184 0.439 0.244 
CPMP P-12 12 0.390 1.053 0.044 0.210 0.538 0.300 
CPMP-envy P=12,d*=0.4km,a=20.0 12 0.399 0.957 0.030 0.173 0.434 0.243 
CPMP P=13 13 0.374 1.053 0.043 0.195 0.521 0.308 
CPMP-envy P=13,d*=0.4km,a=20.0 13 0.386 1.006 0.028 0.168 0.435 0.240 


备注 : CPMP: 容量 约束 P 中 心 问题 ， CPMP-envy: 最 小 妒忌 公平 性 区 位 模型 
2 为 CPMP 和 CPMP-envy 模型 计算 结果 示意 图 。 从 图 中 可 以 看 出 ， 目 标 函 数 增加 最 
小 妒忌 指标 后 ， 设 施 区 位 选择 有 很 大 变化 ，12 个 设施 中 仅 2 个 位 置 相同 。 


图 2 ZY 案例 CPMP( 左 ) 和 CPMP-envy( 右 ) 模 型 计算 结果 示意 图 
Fig. 2 CPMP (left) and CPMP-envy(right) solutions from instance ZY 


案例 ZZ 计算 结果 统计 见 表 4。 表 中 各 栏目 含义 与 表 1 相同 。 对 于 CPMP 模型 (P=64)， 


平均 距离 、 最 大 距离 、 方 差 、 标 准 差 、 变 异 系 数 和 Gini 系数 依次 为 0.727、3.596、0.154、 
0.392、0.540 和 0.293; 增加 最 小 妨 鼠 目标 ， 平 均 距离 增加 了 3.4%， 但 最 大 距离 、 方 差 标准 
差 、 变 异 系 数 和 Gini 系数 减 小 显著 ， 分 别 减少 10.6%、14.3%、7.40%、10.6% 和 6.14%。 设 
施 数 量 减 少 (P=58), CPMP 结果 中 , 除 平 均 距 离 指标 外 , 各 指标 均 变 差 , 增加 最 小 妒忌 目标 ， 
公平 性 指标 明显 改善 ， 设 施 数量 增加 (P=70)， 呈 现 出 类 似 结 

表 4 案例 ZZ 计算 结果 统计 


Tab.4 Computational results from instance ZZ 


模型 。 ”模型 参数 M GNE ME E M A 
CPMP P-58 58 0.758 3.897 0.178 0.422 0.557 0.303 
CPMP-envy P=58,d°=0.75km,a=10.7 58 0.782 3.298 0.147 0.383 0.490 0.271 
CPMP P-64 64 0.727 3.596 0.154 0.392 0.540 0.293 
CPMP-envy P=64,d°=0.75km,a=10.7 64 0.752 3.215 0.132 0.363 0.483 0.275 
CPMP P=70 70 0.692 3.298 0.147 0.383 0.530 0.271 
CPMP-envy P=70,d*=0.75km,a=10.7 70 0.706 2.336 0.105 0.324 0.426 0.258 


备注 : CPMP: 容量 约束 P 中 心 问题 ， CPMP-envy: 最 小 妒忌 公平 性 区 位 模型 
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3 为 CPMP 和 CPMP-envy 模型 P=64) 计 算 结 果 示 意图 。 增 加 最 小 妒忌 目标 后 ， 设 施 
布局 有 明显 的 变化 。CPMP-envy 模型 结果 中 ，36 个 设施 位 置 与 CPMP 模型 结果 相同 。 


图 3 ZZ 案例 CPMP( 左 ) 和 CPMP-envy( 右 ) 模 型 计算 结果 示意 图 
Fig.3 CPMP (left) and CPMP-envy(right) solutions from instance ZZ 


3.3 妒忌 参数 敏感 性 分 析 


本 文 妒忌 指标 计算 与 妒忌 距离 阀 值 参数 必 相 关 。 设 置 这 一 参数 可 能 对 模型 结果 产生 较 大 
的 影响 。 一 般 地 ， 使 用 CPMP 模型 ， 获 得 使 用 服务 的 平均 出 行距 离 ， 参 考 这 一 距离 设置 参数 
d*。 也 可 以 参照 行业 规范 设置 参数 d*, 例如 ，15 分 钟 生活 圈 规 划 推 荐 了 各 类 服务 设施 的 服务 半 
径 ， 可 以 参考 这 一 半径 设置 闽 值 。 本 节 尝 试 设置 不 同 的 妨 忌 距离 冰 值 参数 ， 分 析 这 一 参数 的 
敏感 性 。 针 对 每 个 案例 设置 3 个 参数 ， 进 行 模型 计算 。 计 算 结果 如 表 5 所 示 。 表 5 表明 : (1) 
参数 dq* 对 模型 结果 具有 一 定 的 影响 , 设置 不 同 参数 获得 不 同 的 平均 距离 、 最 大 距离 、 方 差 和 Gini 
系数 ;，(2) 与 原 模型 相 比 ， 增 加 妒忌 指标 后 ， 平 均 距 离 会 有 增加 ， 但 最 大 距离 、 方 差 和 Gini 系数 
指标 均 明 显 改善 ，(3) 需 求 分 布 不 均匀 的 农村 地 区 ， 距 离 变 异 系数 处 于 0.65-0.70 区 间 ， 而 需求 分 
布 较为 均匀 的 城市 地 区 ， 距 离 变 异 系数 处 于 0.43-0.50 区 间 。 

表 5 设置 不 同 妒忌 距离 阀 值 的 CPMP-envy 模型 结果 统计 


Tab. 5 Computational results of minimum envy location models with different envy parameters 


区 域 模型 参数 设施 平均 最 大 方差 标准 差 变异 Gini 


数量 距离 距离 系数 系数 
GY P=22, d*=1.8km,a=4.0 22 2.141 7.002 1.963 1.401 0.654 0.372 
GY P=22, d*=2.0km,a=4.0 22 2.151 8.218 1.988 1.410 0.655 0.370 
GY P=22, d*=2.2km,a=4.0 22 2.102 8.218 2.109 1.452 0.691 0.391 
ZY P=11, d*-0.35km,a-20.0 11 0.422 1.006 0.035 0.187 0.443 0.248 
ZY P=11, d*=0.40km,a=20.0 11 0.420 1.006 0.034 0.184 0.439 0.244 
ZY P=11, d*=0.45km,a=20.0 11 0.424 1.006 0.034 0.184 0.435 0.242 
ZZ P=64, d*=0.70km,a=10.0 64 0.742 3.896 0.136 0.369 0.497 0.273 
ZZ P=64, d*=0.75km,a=10.0 64 0.752 3.215 0.137 0.370 0.492 0.275 
ZZ P=64, d*=0.80km,a=10.0 64 0.754 2.499 0.128 0.358 0.474 0.266 


3.4 与 其 他 公平 性 模型 比较 


考虑 到 距离 方差 被 推荐 为 公平 性 区 位 问题 的 常用 指标 5149， 本 证 尝试 使 用 MVLP 和 
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MVLP-MD 模型 对 案例 GY 进行 求解 。.MVLP 和 MVLP-MD 模型 使 用 CPLEX 12.6 优化 器 直 
接 求解 。 表 6 为 MVLP 和 MVLP-MD 模型 部 分 计算 结果 。 可 以 看 出 ，MVLP 模型 结果 中 ， 
距离 方差 和 Gini 系数 很 小 ， 但 平均 距离 很 大 ; MVLP-MD 能 够 根据 最 大 均 距 约束 获得 不 同 
的 设施 布局 ， 随 参数 qd’ 的 增 大 ， 方 差 和 Gini 系数 逐步 减 小 。 

表 6 案例 GY 计算 结果 统计 


Tab. 6 Computational results from instance GY 


模型 模型 参数 2" D a WM m 
MVLP P=23 23 16.999 18.967 0.159 0.399 0.023 0.013 
MVLP-MD P=23,d'=8.0km 23 8.000 10.157 0.206 0.454 0.057 0.031 
MVLP-MD P=23,d'=3.0km 23 3.000 7.163 0.816 0.903 0.301 0.163 
MVLP-MD P=23,d'=2.0km 23 2.000 8.218 1.976 1.406 0.703 0.396 


备注 : MVLP: 方差 最 小 区 位 模型 ，MVLP-MD: 最 大 距离 约束 的 方差 最 小 区 位 模型 


值得 注意 的 是 ，MVLP 模型 计算 复杂 度 极 高 ， 案 例 GY 在 24h 计算 时 间 内 ， 仍 不 能 获得 
最 优 解 ， 且 方差 目标 严重 扭曲 效率 指标 。 图 4 左 图 为 MVLP 模型 结果 ， 绿 色 方 块 表 示 设 施 
位 置 ， 设 施 服务 区 (设施 服务 需求 点 集合 ) 使 用 不 同 颜色 表 示 。 模 型 结果 具有 以 下 特征 : 设施 
服务 区 凌乱 ; 设施 分 布 于 需求 密度 稀 疏 区 域 ; 设施 并 不 服务 临近 需求 点 , 而 是 服务 距离 16.0- 
18.0km 范围 内 的 部 分 需求 点 。 虽 然 MVLP 模型 结果 明显 不 符合 常理 ， 但 公平 性 指标 很 好 ， 
方差 和 Gini 系数 分 别 为 0.159 和 0.013. 


4 GY 案例 MVLP( 左 ) 和 MVLP-MD( 右 ) 模 型 计算 结果 示意 图 
Fig. 4 MVLP (left) and MVLP-MD (right) solutions from instance GY 
图 4 右 图 为 MVLP-MD 模型 (d'=3.0km) 结 果 。 与 MVLP 相 比 ，MVLP-MD 模型 计算 复 


杂 度 大 幅 降 低 ， 倾 向 于 在 需求 密度 较 高 区 域 选 择 设施 区 位 ; 但 在 部 分 区 域 ， 供 需 空 间 关 系 仍 
然 被 扭曲 。 当 参数 qd' 接 近 CPMP 模型 最 优 解 时 (最 优 距 离 为 1.948km)，MVLP-MD 与 CPMP 
结果 接近 。 

图 5 为 距离 参数 d'( 横 轴 ， 单 位 km) 与 最 小 方差 ( 纵 轴 ) 关 系 。 当 qd' <1.948km 时 ，MVLP- 
MD 模型 无 可 行 解 ， 增加 d' 数 值 ， 平 均 距 离 略 小 于 必 ， 但 最 小 方差 显著 下 降 ， 当 d' >3.0km 
时 , 最 小 方差 下 降幅 度 趋 缓 ; 4d > 17.0km 之 后 , 最 小 方差 走势 改变 , 开始 呈现 增加 趋势 ; 
当 d' > 22.0km 之 后 ， 最 小 方差 呈现 显著 增加 趋势 。GY 案例 模型 计算 表明 : MVLP-MD F 
均 距 离 与 最 小 方差 呈现 出 长 平底 U 型 曲线 。 
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图 5 GY 案例 MVLP-MD 模型 平均 距离 ( 横 轴 ) 与 最 小 方差 ( 纵 轴 ) 之 关系 


Fig. 5 Travel distance mean (horizontal) versus variance (vertical) of instance GY 


空间 可 达 性 分 析 是 评价 和 改善 服务 设施 服务 质量 及 其 公平 性 的 重要 手段 。 基 于 2SFCA 
可 达 性 方法 ，2SO4SAI 越 来 越 受 到 关注 。 本 节 采 用 该 方法 进行 案例 测试 ， 首 先 使 用 PMP E 
型 选择 设施 区 位 , 再 参考 文献 2 构建 最 小 化 2SFCA 指标 方差 二 次 模型 进行 公平 性 资源 分 配 。 
2SO4SAI 方法 求解 过 程 如 下 : 首先 使 用 CPLEX 12.6 求解 PMP 模型 ， 再 参考 文献 5 建立 二 
次 规划 模型 ， 使 用 CPLEX 12.6 求解 。 

对 于 县 级 市 案例 GY, 选择 设施 数量 P=23, 求解 PMP 获得 设施 区 位 及 其 服务 区 (设施 服 
务 需求 点 的 集合 )。 使 用 服务 出 行距 离 、 最 大 距离 、 方 差 和 Gini 系数 分 别 为 1.866、8.218、 
1.898 和 0.414。 设置 资源 总 量 为 900000, 设施 服务 距离 闵 值 为 4.0km 和 5.0km。 资 源 分 配 结 
果 如 图 6 上 左 所 示 。 图 中 ， 设 施 位 置 的 三 个 直方 图 自 左 向 右 依次 为 设施 服务 区 总 需求 量 、 
4.0km 距离 闷 值 资源 分 配 量 和 5.0km 距离 闵 值 资源 分 配 量 。 可 以 看 出 : 资源 分 配 结果 受 距离 
闵 值 参数 影响 显著 ; 资源 需求 量 与 分 配 量 有 比较 明显 的 差异 ; AA E. 外围 区 域 的 设施 资源 
分 配 相 对 于 需求 偏 少 ， 而 中 部 区 域 ， 资 源 分 配 相 对 于 需求 偏 高 。 

对 于 城市 案例 ZY， 选 择 设施 数量 P=12， 求 解 PMP 获得 设施 区 位 及 其 服务 区 。 使 用 服 
务 出 行距 离 、 最 大 距离 、 方 差 和 Gini 系数 分 别 为 0.388、1.053、0.042 和 0.295。 设 置 资源 总 
量 为 4000， 设 施 服务 距离 阔 值 为 0.6km、0.8km 和 1.0km。 资 源 分 配 结果 如 图 6 上 右 所 示 。 
与 案例 GY 相似 , 资源 分 配 结果 受 距 离 闵 值 参数 影响 显著 ; 资源 需求 量 与 分 配 量 有 较为 明显 
的 差异 ， 但 相对 差异 小 于 案例 GY。 

对 于 大 型 城市 案例 ZZ， 选 择 设施 数量 P=60， 求 解 PMP 获得 设施 区 位 及 其 服务 区 。 使 
用 服务 出 行距 离 、 最 大 距离 、 方 差 和 Gini 系数 分 别 为 0.738、3.897、0.143 和 0.279。 设 置 资 
源 总 量 为 4500000, 设施 服务 距离 阔 值 为 1.2km、1.5km 和 1.8km. 资源 分 配 结果 如 图 6 下 所 
示 。 图 中 , 设施 位 置 的 四 个 直方 图 最 左 向 右 依次 为 设施 服务 区 总 需求 量 、 三 个 距离 闹 值 资源 
分 配 量 。 可 以 看 出 : 资源 分 配 结果 受 距 离 闵 值 参 数 影响 显著 ; 资源 需求 量 与 分 配 量 有 比较 明 
显 的 差异 ; 外 围 区 域 的 设施 资源 分 配 相 对 于 需求 偏 少 , 而 中 部 区 域 , 资源 分 配 通常 高 于 需求 。 
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图 6 案例 GY( 上 左 )、ZY( 上 右 ) 和 ZZ( 下 )2SO4SAI 模型 计算 结果 示意 图 
Fig. 6 2SO4SAI solutions from instances GY (upper left), ZY (upper right) and ZZ (bottom) 


结合 三 个 案例 的 需求 分 布 特征 , 可 以 发 现 : 2SO4SAI 受 可 达 性 模型 距离 参数 影响 显著 ; 


在 人 口 密度 高 的 地 区 ， 相对 于 需求 量 , 模型 倾向 于 为 设施 分 配 更 多 的 资源 ; 而 在 需求 分 布 稀 
MAER, 相对 于 需求 量 , 模型 倾向 于 为 设施 分 配 更 少 的 资源 ,作者 认为 , 有 必要 针对 2SO4SAI 
方法 设计 若干 测试 案例 ， 系 统 地 进行 参数 敏感 性 分 析 、 供 需 关 系 分 析 、 公 平 性 指标 分 析 ， 进 
一 步 验 证 该 方法 的 基本 特征 及 其 可 用 性 。 


4 讨论 与 结论 T6 


H 1990s 初 ,学 者 逐步 使 用 公平 性 指标 构建 公平 性 区 位 模型 。 第 用 公平 性 指标 包括 距离 
方差 、 基 尼 系 数 、 期 望 值 偏差 、 可 达 性 方差 等 。 然 而 ， 将 这 些 指标 用 于 设施 选 址 ， 模 型 优化 
结果 差异 巨大 。 除 2SOASAT 有 少量 案例 应 用 研究 外 ， 鲜 见 最 小 化 距离 方差 、 基 尼 系 数目 标 
函数 的 实际 应 用 。 出 现 这 一 情况 的 原因 在 于 : 距离 方差 、 基 尼 系 数目 标 会 严重 扭曲 空间 可 及 
性 指标 ， 导 致 模型 结果 不 可 用 。 因 此 ， 距 离 方 差 、 基 尼 系 数 可 以 用 于 评价 设施 布局 规划 的 空 
间 公 平 性 , 但 不 宜 作为 优化 目标 。 对 于 离散 设施 选 址 ， 本文 MVLP-MD 模型 实验 表明 : 设施 
数量 保持 不 变 ， 出 行距 离 均值 与 最 小 距离 方差 呈现 出 一 条 长 平底 U 型 曲线 (图 3)。 随 出 行距 
离 均 值 增加 ， 距 离 方差 快速 下 降 ， 但 下 降 趋 势 逐 步 变 组 ， 出行 距离 均值 达到 某 一 阐 值 后 ， 方 
兰 处 于 很 低 水 平 ; 出 行距 离 均 值 达到 茶 一 很 大 的 病 值 时 , 方差 又 开始 增加 。 在 曲线 底部 的 距 
离 区 间 内 ， 方 差 指标 表现 优异 ， 但 在 这 一 区 间 ， 出 行距 离 被 严重 扭曲 。 
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基于 图 5 所 示 的 出 行距 离 与 距离 最 小 方差 之 间 的 关系 ， 有 三 个 重要 的 启示 : (1) 区 位 模 
型 不 能 将 距离 方差 作为 优化 目标 , 因为 U 型 曲线 方差 较 小 区 域 , 出 行距 离 已 经 被 严重 扭曲 ; 
(2) 距 离 方差 目标 配合 出 行距 离 约束 仍旧 不 具有 可 行 性 , 出 行距 离 设 置 较 大 时 , 模型 结果 无 实 
际 意义 ， 即使 出 行距 离 设置 较 小 (U 型 曲线 最 左 的 区 间 ), 部 分 区 域 的 出 行距 离 被 会 被 扭曲 (如 
图 1 右 图 所 示 ); (3) 距 离 方 差 也 不 宜 作 为 多 目标 优化 的 目标 沙 数 ， 原 因 在 于 模型 Pareto 解 集 
中 ， 大 多 数 解 的 效率 指标 被 严重 扭曲 。 

基于 2SFCA 空间 可 达 性 模型 进行 公平 性 资源 分 配 ， 三 个 区 域 案 例 的 实验 结果 表明 : 多 
数 设施 服务 区 内 会 出 现 供需 失衡 。 作 者 认为 ， 出 现 该 现象 的 根本 原因 在 于 2SFCA 可 达 性 模 
型 本 身 的 局 限 : 可 达 性 指标 对 于 距离 衰减 模型 选择 、 模 型 参数 设置 敏感 ， 计 算 结果 的 物理 意 
义 不 够 清晰 中。 本 质 上 ， 一 个 服务 系统 的 供需 关系 很 难 使 用 一 个 距离 衰减 模型 进行 准确 表 
达 ， 因 服务 供需 关系 不 仅 受到 距离 影响 ， 也 与 需求 空间 分 布 特征 、 供 应 分 布 特征 密切 相关 。 

本 文 将 最 小 方差 指标 修订 为 最 小 “妒忌” 指标 。 该 修订 基于 一 个 假设 : 使 用 服务 的 出 行距 
离 较 近 时 ， 需 求 者 不 会 有 不 公平 感知 ， 或 者 说 不 会 嫉妒 其 他 需求 者 ;使 用 服务 出 行 较 远 而 不 
够 便捷 时 ， 需 求 者 会 产生 不 公平 感知 ， 或 者 说 会 嫉妒 其 他 需求 者 ， 不 公平 测度 或 者 妒忌 测度 
与 出 行距 离 相关 。 基 于 这 一 假设 ， 本 文 嫉妒 指标 仅 计 算 距 离 超 过 某 一 阔 值 的 需求 者 。 这 一 指 
标 与 经 典 模型 的 成 本 指标 不 冲突 , 乘 以 一 个 特定 系数 后 , 妒忌 指标 量 纲 与 出 行 成 本 量 纲 也 能 
够 相 一 致 ， 从 而 易于 集成 在 经 典 模型 中 。 本 文 三 个 区 域 案 例 的 实验 表明 : 最 小 妒忌 指标 能 够 
用 于 改进 CPMP 模型 的 目标 函数 ， 目 标 函 数 改 进 后 ， 出 行距 离 略 有 增加 (增幅 约 2~7%)， 同 
时 ， 最 大 距离 、 方 差 、 标 准 差 、 变 异 系数 和 Gini 系数 等 公平 性 指标 显著 改善 ， 当 设施 数量 
较 少 时 ， 公 平 性 指标 改善 尤为 显著 。 

作者 认为 , 最 小 妒忌 指标 在 公平 性 区 位 问题 建 模 中 具有 理论 价值 , 在 公共 服务 设施 布局 
规划 中 具有 应 用 潜力 。 进 一 步 的 研究 工作 包括 : 从 理论 上 验证 公共 服务 妒忌 假设 能 够 成 立 ， 
并 建立 妨 忌 度量 模型 ， 推 荐 合适 的 模型 结构 及 参数 取 值 ， 进一步 改进 模型 求解 算法 ， 能够 高 
效 高 质量 地 求解 新 模型 ; 选择 更 多 的 实际 案例 测试 并 验证 最 小 妒忌 指标 可 用 性 ; 探索 公共 服 
务 最 小 妒忌 指标 、 供 给 成 本 指标 、 空 间 可 及 性 指标 、 公 平 性 指标 之 间 可 能 存在 的 规律 。 
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A minimum-envy location problem for spatial equality of public service: 


model formulation and case analysis 


Yunfeng Kong!”, Shiyan Zhai!” 

(1 Key Laboratory of Geospatial Technology for the Middle and Lower Yellow River Regions, Ministry of 
Education, Henan University, Kaifeng, China 47500; 2 College of Geography and Environmental Science, Henan 
University, Kaifeng, China 47500) 

Abstract: Location problems have been widely applied to service planning of public health, 
compulsory education, emergency management, and delivery logistics. However, the mainstream 
location models are usually to optimize the efficiency objectives such as travel cost, facility cost 
and the number of customs to be served, rather than the equality objectives. A few location models 
aim to optimize one of the equality measures, such as the variance of distances, the deviation of 
distances, the Gini coefficient between the travel distances, and the variance of spatial accessibility 
indexes. However, the facility locations, capacities and their service areas can be easily distorted by 


most equality-oriented objective functions. In this paper, a spatially envy objective function for 
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service equality is proposed to overcome the shortcomings of commonly used equality functions. 
The envy value of customers at a location is determined by their travel distance that beyond a 
predefined distance. The envy function can be added to mainstream location models in a weighted 
manner. As a result, the capacitated p-median problem (CPMP) is enhanced as CPMP-envy. The 
original and improved models were tested on three large instances. Case experiments show that the 
equality measures, such as maximum travel distance, variance of distances, coefficient of variation, 
and Gini coefficient between travel distances, can be substantially improved by minimizing the 
weighted sum of spatial envy and travel cost. It is argued that the envy indicator has theoretical and 
practical potentials in facility planning towards spatial equality of public service. 

Keywords: public service; location model; spatial equality; minimum variance; minimum envy; 


case analysis. 
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